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The final state of Pb + Pb reactions at CERN SPS has been reconstructed with the 
Buda-Lund hydro model, by performing a simultaneous fit to NA49, NA44 and WA98 
data on particle correlations and spectra. 

1. Introduction 

In refs. [0,0 it was observed, for the first time, that the parameters of particle emitting 
sources can be determined only if a simultaneous analysis of the momentum distributions 
and two-particle correlation functions is performed. Simultaneous fitting of particle cor- 
relations and spectra were reported in refs. 0-0]. Here, we determine the reconstructed 
space-time picture of particle emission in Pb + Pb collisions at CERN SPS by fitting 
simultaneously the NA44 and NA49 published data [|,[| on two-particle correlations and 
single-particle spectra at Pb + Pb 158 AGeV central reactions at CERN SPS. Preliminary 
data [jTOj from WA98 experiment are also used to check the reliability and the consistency 



of the fit results. 



2. Buda-Lund hydrodynamic model 

The Buda-Lund hydro parameterization characterizes with means and variances the 
local temperature, flow and chemical potential distributions of a cylindrically symmet- 
ric, finite hydrodynamically expanding system. The four- velocity u^(x) of the expanding 
matter is given by a scaling longitudinal Bjorken flow appended with a linear transverse 
flow, characterized by its mean value, (u t ). A Gaussian shape of the local density distri- 
bution is assumed both in the transverse plane and in space-time rapidity. The changes 
of the inverse temperature are characterized with means and variances. The freeze-out 
hypersurface is characterized by a mean freeze-out (proper) time tq and a duration param- 
eter At, the variance of the freeze-out propertime distribution. The following emission 
function S c (x,p) applies: 
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exp l^r-^)j +s 

p^dH^ix) = m t cosh[r/ — y]H(r)dTTodr] dr x dr y , (2) 
u^(x) = (cosh [77] cosh [rjt], sinh[77 t ] r x /r t , sinh[^ t ]r y /r t , sinh [77] cosh [r) t ]) , (3) 

Ms) ^ r * + r l (V- Vo) 2 T -T r = AT To-T t _ AT 
T(x) " T 2i?^ 2Ar/ 2 ' T r T T t T 



where is the chemical potential and T(x) is the local temperature, the subscript c 
refers to the core of collision (surrounded by a halo of long lived resonances), the transverse 



mass is m t = ^m? + p x + p y , the rapidity y and the space-time rapidity 77 are defined 
as y = 0.5 log [(E + p z )/(E — p z )] and r\ = 0.5 log [(t + z)/(t — z)], Rg stands for the 



transverse geometrical radius of the source, sinh [77^ = (u t )-^- with r t = \/r%+r% . The 

Rg * 

central temperature at mean freeze-out time is denoted by To = T{r x = r y = 0; r = To) ■ 
With the surface temperature T r = T{r x = r y = Rg, t = r ) and the temperature after 
freeze-out, T t = T(r x = r y = 0; r = tq + a/2At), the relative transverse and temporal 
temperature decrease are introduced, see refs. ||,|2|,[l| for further details. 

2.1. Single particle spectra and two particle correlations 

The total invariant single particle spectrum in y rapidity and transverse mass mt is 
determined analytically using a saddle-point approximation, 

3 EVC , (5) 



2irm t dm t dy (27r) / u^(x)p 

exp 1 




Tix) 

where E stands for the average energy, V for the average volume of the effective source 
of particles with a given momentum p. The correction factor C includes the effective 
intercept parameter \*(y,mt) of the two particle (or Bose-Einstein) correlation function, 
that controls the core ratio in the particle production in the core/halo picture. In Ref. [TJ, a 
Boltzmann approximation (s = 0) to the above Invariant Momentum Distribution (IMD) 
as well as to the Bose-Einstein correlation function (BECF) was obtained in an analytic 
way. The analytical formulas for the BECF and IMD, as were used in the fits, have been 
summarized in refs. JI|,0,||,Q. Bose-Einstein or Fermi-Dirac statistics (s = ±1) was used 
in the analytic expressions fitted to single particle spectra. 



3. Fitting NA49, NA44 and WA98 Pb + Pb data 

The kinematic parameters of the Buda-Lund model are fitted simultaneously to IMD 
and HBT radii measured by the CERN NA49, NA44 and WA98 experiments in central 
Pb + Pb collisions at 158 AGeV. Core/halo correction oc is applied and the corre- 

sponding errors are propagated properly. Due to the these conditions unique minima are 
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found, a good x 2 /NDF is obtained for all reactions, and the strongly coupled, normal- 
ization sensitive {^f-)t an d Ar parameters are determined. On Figure 1 and 2, the fits to 
measured data are shown together with the published data. Note that the first 5 points 
of the NA44 pion spectrum were contaminated 11 and were not included in the fit. 
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Figure 1. Simultaneous fits to NA49 parti- Figure 2. Simultaneous fits to NA44 parti- 
cle spectra and HBT radius parameters. cle spectra and HBT radius parameters. 



Fits to preliminary data of WA98 experiment provide source parameter values and 
errors similar to those obtained by NA49 and NA44. The normalizations of the NA44 
pion spectrum and the WA98 h~ spectrum had to be fixed manually to that of NA49. 

The hypothesis that pions, kaons and protons are emitted from the same hydrody- 
namical source is in a good agreement with all the fitted data. The fit parameters are 
summarized in Table 1, shown with statistical errors, only. The reconsructed space-time 
emission function S(x) (which is the source function S(x,p) integrated over the momen- 
tum p) is shown on Figure 3 and 4. 



4. Conclusions 

We find that the NA49, NA44 and WA98 data on single particle spectra of h~ , identified 
7r, K and p as well as detailed rapidity and m t dependent HBT radius parameters are 
consistent with each other. The final state of central Pb + Pb collisions at CERN SPS 
corresponds to a cylindrically symmetric, large (Rq = 7.1 ± 0.2 fm) and homogenous 
(To = 139 ± 6 MeV) fireball, expanding three-dimensionally with (u t ) = 0.55 ± 0.06. A 
large mean freeze-out time, r = 5.9 ± 0.6 is found with a short duration of emission. 
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Table 1 

Source parameters from simultaneous fitting of NA49, NA44 and preliminary WA98 par- 
ticle spectra and HBT radius parameters with the Buda - Lund hydrodynamical model. 





NA49 


NA44 


WA98 


Averaf 


led 


Parameter 


Value 


Error 


Value 


Error 


Value 


Error 


Value 


Error 


To [MeV] 


134 


± 3 


145 


± 3 


139 


± 5 


139 


± 6 


(ut) 


0.61 


± 0.05 


0.57 


± 0.12 


0.50 


± 0.09 


0.55 


± 0.06 


Rg [H 


7.3 


± 0.3 


6.9 


± 1.1 


6.9 


± 0.4 


7.1 


± 0.2 


t [fm/c] 


6.1 


± 0.2 


6.1 


± 0.9 


5.2 


± 0.3 


5.9 


± 0.6 


At [fm/c] 


2.8 


± 0.4 


0.01 


± 2.2 


2.0 


± 1.9 


1.6 


± 1.5 


Arj 


2.1 


± 0.2 


2.4 


± 1.6 


1.7 


± 0.1 


2.1 


± 0.4 


(f)r 


0.07 


± 0.02 


0.08 


± 0.08 


0.01 


± 0.02 


0.06 


± 0.05 


0.16 


± 0.05 


0.87 


± 0.72 


0.74 


± 0.08 


0.59 


± 0.38 


X 2 /NDF 


163/98 


= 1.66 


63/71 


= 0.89 


115/108 


, = 1.06 


1.20 



Pb + Pb at CERN SPS 



Pb + Pb at CERN SPS 



S(z,t) 





S(x,y) 



Figure 3. The reconstructed source function Figure 4. The reconstructed source function 
S(t, z, x — 0, y — 0) in the (t, z) plane. S(x, y) at the mean freeze-out time. 
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